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A laterally confined thin elastic sheet lying on a liquid substrate displays regular undulations, 
called wrinkles, characterized by a spatially extended energy distribution and a well-defined wave¬ 
length A. As the confinement increases, the deformation energy is progressively localized into a 
single narrow fold. An exact solution for the deformation of an infinite sheet was previously found, 
indicating that wrinkles in an infinite sheet are unstable against localization for arbitrarily small 
confinement. We present an extension of the theory to sheets of finite length L, accounting for the 
experimentally observed wrinkle-to-fold transition. We derive an exact solution for the periodic de¬ 
formation in the wrinkled state, and an approximate solution for the localized, folded state. We find 
that a second-order transition between these two states occurs at a critical confinement Af = A^/L. 

PACS numbers: 46.32.-|-x, 46.70.-p, 46.70.De, 


I. INTRODUCTION 

Morphological transitions are often induced by con¬ 
finement or by spatially constrained growth. Structures 
emerge spontaneously when the energy injected into a 
system through the confinement process ceases to dis¬ 
tribute uniformly. These phenomena are observed in var¬ 
ious contexts ranging from the folding of geological lay¬ 
ers [1] to patterns in biological membranes and monolay¬ 
ers [2-7] and the formation of fingerprints [8, 9]. Besides 
the initial morphology, occurring for small confinement, 
various subsequent transitions may be observed as the 
confinement increases. 

Many model systems have been developed to study 
the influence of confinement on morphological transi¬ 
tions [10-16] using rods [17-19] or sheets resting on some 
substrate [20-31]. Among them, the experimental model 
system of a thin elastic sheet, lying on a fluid sub¬ 
strate [32-41] and subjected to in-plane uniaxial com¬ 
pression, deserves special interest [42-49] . Beyond a cer¬ 
tain critical confinement. Aw, which vanishes for an in¬ 
compressible sheet, the sheet buckles, displaying regular 
undulations with a characteristic wavelength A over its 
entire length [2, 42] (Fig. 1). Beyond a certain confine¬ 
ment, Ap, another transition occurs, where the wrinkles 
start attenuating except near the center of the sheet, 
and the deformation energy gets localized into a single 
fold [42]. 

In the case of an infinitely long sheet, an exact solution 
describing the complete evolution of the morphology with 
increasing confinement proves that such a wrinkle-to-fold 
transition, strictly, should never occur [45], i.e. Ap —>■ 
Aw- The infinite-sheet morphology is always localized. 
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FIG. 1: (Color online) Schematic morphological evolution of 
a compressed finite sheet floating on a liquid. The scheme 
shows the emergence of a wavelength A upon small confine¬ 
ment, A > Aw, the growth in amplitude of the wrinkle state 
for intermediate compression, and the transition to a fold 
state beyond some critical confinement, Ap. Here t and ft 
are the tangent and normal to the sheet surface, respectively. 
(j>{s) is the angle between the local tangent and the horizontal 
direction x, and the arclength along the sheet is parametrized 
by s. h{s) is the vertical elevation of the sheet. 


but the localization decay length diverges as the buckling 
threshold Aw is approached. 

Even if an infinitely long sheet is a useful idealization 
of real systems allowing a good agreement with experi¬ 
ment [46] , the apparent discrepancy between theory and 
experiment concerning the existence of a wrinkle-to-fold 
transition should be resolved. Experiments are obviously 
performed with finite sheets. We demonstrate that a 
second-order wrinkle-to-fold transition does occur for fi¬ 
nite sheets. In Sec. Ill we construct an exact periodic 
solution for wrinkles of arbitrary amplitude in a finite 
sheet. An exact solution for a localized deformation in 
a finite sheet has evaded us. Hence, in Sec. IV, we de¬ 
rive an approximate localized solution using a multiple- 
scale analysis. In Sec. V, we identify the order parame¬ 
ter of the transition in the context of the Landau theory 
of second-order transitions. We show how both periodic 
and localized solutions match at the critical confinement. 
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Ap = /L. In Sec. VI, we discuss the experimental im¬ 

plications of our theory, its limitations, and future exten¬ 
sions. 


II. SYSTEM AND GOVERNING EQUATION 


The system studied here is composed of a thin incom¬ 
pressible elastic sheet (Aw = 0) of length L, width W, 
and bending modulus B, lying on a fluid of mass density 
p. The sheet is uniaxially confined by a distance A along 
the x-axis and deforms in the xy plane. The shape of the 
sheet is described by the parametric equation 


x(s) = / cos(j){s)ds', 

J-L/2 

h{s) = / siii4>{s')ds', 

J-L/2 


(la) 

(lb) 


where 4>(s) is the angle between the local tangent to the 
sheet and the x-axis at a given arclength s, see Fig. 1. 
The total energy of the system, E, is composed of the 
bending energy of the sheet, = {WBj2) 
and the deformation energy of the substrate, Eg = 
{Wpg/2) cos (j)ds, where the dot denotes an s 

derivative [42, 45]. The displacement along the direction 
of confinement is given by 


A = 


fL/2 

/ (1 —C0S(^)ds, 

J-L/2 


( 2 ) 


and is related to the applied load necessary to confine the 
sheet by P = dE/dA. In the following, except where it is 
explicitly mentioned, we use units such that the energy is 
rescaled by B, and lengths are rescaled by (Bjpg)^!^ = 
A/27r. As a result, the applied load P is rescaled by 
{Bpgf!'^. 

To find the equilibrium shape, one should minimize the 
total energy under appropriate constraints. This can be 
reformulated as a dynamical problem [45] with an action 
S = ^)’ ■''^here 


£ = + coscj)—P{1—cos (p—A/L) — Q (s) {sin (j)—h). 

( 3 ) 

In Eq. (3), P and Q{s) are Lagrange multipliers in¬ 
troduced to take into account, respectively, the global 
constraint (2) and the local geometrical constraint be¬ 
tween h and (p. The conjugate momenta are defined as 
= dCpdcp = (p and pu = dCjdh = Q and are used 
to construct the Hamiltonian % = p^cp + phh — C. Since 
C has no explicit dependence on s, the Hamiltonian is a 
constant for a given displacement A, 

Ti. = ^p^+Ph sin((i— cos(()-|-P(l —cos(() —A/L) = C. 

( 4 ) 


Hamilton’s equation, p^p = —dH/dtp, yields the following 
equation: 


0-1-(/i^/2-I-P)sin(()-|-p?iCOS(() = 0. (5) 

Eliminating ph between Eqs. (4) and (5), we obtain 

0sin0 — ^0^ COS0 -|--I-P — Pcos 0 = 0, (6) 

where P = P{\ — A/L) — C is a shifted value of the 
load, dependent on boundary conditions. Differentiat¬ 
ing Eq. (6) once, we obtain the well known equation for 
Euler’s elastica, 

'4> + 4> + h = Q, (7) 

where the local normal force exerted on the elastica is 
given by the hydrostatic term h. A second differentiation 
gives the equation governing the system evolution, 

3 . .. 

0 -I--0^0-I-P0-I-sin0 = 0. (8) 

Finally, Eq. (6) is used to relate the physical load P to 
the shifted value P. For hinged boundary conditions, 
where h and h (and also 0) vanish at the boundaries, we 
have 


P = Pcos0(±L/2) - 0(±L/2)sin0(±P/2). (9) 

As a result, once the solution 0(s) is obtained, the phys¬ 
ical load P can be computed from P. Note that in the 
limit of an infinite sheet, where 0(±P/2) —)■ 0, P and P 
coincide. 

The total rescaled energy per unit length of the system 
reads 


P[0(s)] — ^ J ■ (10) 

The energy, the displacement A and the equation giving 
the equilibrium shape of the sheet can also be written in 
terms of h instead of 0 using h = sin 0: 


E[h{s)] 

A[/i(s)] 


f.L/2 


ds 




J-L/2 

L/2 


1 - yp 




ds 1 — y 1 — , 


h2^,(ll) 

( 12 ) 


J-L/2 \ / 

The shape of the sheet is determined from minimization 
of either P, given the displacement A, or the function. 


given the load P. 


G = E-PA, 


(13) 
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III. EXACT PERIODIC SOLUTIONS 
A. General results 

We first study general solutions of Eq. (8) without 
specifying the boundary conditions. In order to construct 
a periodic solution of Eq. (8), we recall the connection 
existing between this equation and the dynamics of a 
physical pendulum [45, 52]. For this purpose, we con¬ 
sider the total energy Up of a pendulum, 

C7p = ^-t g^(l - COS0), (14) 

where q is the natural frequency of the pendulum, and k 
is a constant determined by the pendulum’s total energy 
or, equivalently, by boundary conditions. (In the anal¬ 
ogous elastic sheet these two parameters are related to 
the natural wrinkling wavenumber, 27r/A, and the total 
displacement A.) The equation of motion is obtained by 
taking the variation of Up in the above equation, 

(j) + smcj) = 0. (15) 

Differentiating twice this last equation and eliminating 
cos ^ and sin (/) using respectively Eqs. (14) and (15), 
we obtain 

^ + (l - 2A:^) g^((i = 0. (16) 

Adding Eq. (15) multiplied by q~^, we finally obtain the 
equation describing the morphology of a confined floating 
sheet, see Eq. (8): 


3 . .. 

-I- P(j) -1- sine/ = 0, 

(17a) 

g2(l_2fc2)+g-2. 

(17b) 


where £{x,k) is the Jacobi epsilon function [50, p.562]. 
The quantities sq and g are fixed to satisfy the bound¬ 
ary conditions whereas the parameter k is related to the 
confinement ratio, A/L. 


B. Hinged sheets 

For hinged boundary conditions both h and h vanish 
at s = Ah 12. The second derivative of the sheet profile, 
Eq. (19), is given by 

h{s) = 2kq cn(g(s -|- Sq), k) (2 dn^(g(s -|- sq)) ~ l) • 

( 21 ) 

The dn function varies between (1 — and 1; conse- 

quently, so long as k < ^.1^/2, the last factor in Eq. (21) 
does not vanish. As shown below, the relevant values of 
k are smaller than 1/3 provided L > 3X. Larger values of 
k lead to periodic solutions unstable against localization. 

Due to the periodicity of the solution, there exists an 
infinite number of possible solutions depending on the 
number of nodes. From Eq. (19), h(L/2) = 0 provided 
that gL/2-|-gso = —qL/2+qso+4:niK{k). FromEq. (21), 
’h{L/2) = 0 provided that gL/2 -|- gso = (2n2 -I- VjK{k), 
with Til and n 2 positive integers. Due to the definite 
parity of the solution (either symmetric or antisymmet¬ 
ric), the condition h{—L/2) = 0 is then automatically 
satisfied. We have the following two possibilities, (a) 
If So = 0, then h(L/2) = 0 automatically, for any g, 
due to the even parity of the cn function in Eq. (19). 
We are left with the condition for h{L/2) = 0 which 
gives qL = 2{2n2 + l)K{k). (b) If sq ^ 0, then the two 
conditions above must be satisfied simultaneously, giving 
qL = AniK{k) and gso = [2(n2 — ni) + l]K{k). Thus, 
combining these two results, we have 


Therefore, any solution of Eq. (15) is also a solution of 
Eq. (8) with P given by Eq. (17b). 

It is well known that Eq. (15) admits a periodic so¬ 
lution for 0 < fc < 1 in terms of Jacobian elliptic func¬ 
tions [50, p.549], [56]: 

(j){s) = 2 arcsin [k sn(g(s -|- sq), k)]. (18) 


L 


where p is a strictly positive integer. Finally, the solution 
reads 


The profile h{s) = sin ^(s')ds' is thus given by 

2k 

h{s) = — [cn(g(s -|- Sq), k) - cn(g(-L/2 -|- sq), k)] . 

q 

(19) 

This solution has a periodicity of AK{k) /g, where K(x) is 
the complete elliptic integral of the first kind [50, p.487]. 
Finally, the horizontal displacement is given by 

fL/2 

A = 2k^ / sn^{q{s + So), k) ds 
J-L/2 

= 2L-^{£:(g(s + So),fc)}^:i%, (20) 


(j){s) = 2 arcsin [ksn{qs + eK{k), k )], 
2k 

h{s) = —cn(gs -|- eK{k), k), 

q 

where e = {1 + (—1)”)/2 and 


2nK{k) 

z 


n= 1,2,3... 


(24a) 

(24b) 


(25) 


Symmetric solutions correspond to n of odd parity and 
antisymmetric solutions are obtained with an even parity 
of n. The number of nodes of the solution is equal to n+1 
(counting the two nodes at the boundaries). 
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FIG. 2: (Color online) Evolution of the exact periodic solu¬ 
tion for different values of A/A and L/A = 3.5 {N = 7), where 
x{s), h{s) are given by Eqs. (la) and (24b) respectively and 
(j>{s) is given by Eq. (24a). The corresponding evolution for 
the numerical solution of Eq. (8) is also shown. P and Pp 
correspond to the numerical and the periodic solutions re¬ 
spectively (P = Pp for A/A < 0.2825). The wrinkle-to-fold 
transition occurs at A/A ~ 0.2825. 



C. Pressure, displacement, and amplitude 


From Eq. (24a), we have 

cos(/)(L/2) = 1 - 2fc^ ^(L/2) = 

sin(/i(L/2) = ±2fc-\/l — (26) 

The expression for the applied load is then obtained from 
Eq. (9) together with Eq. (17b): 

P ^ q'^ + {l-2k^)q-^. (27) 


From this infinity of possible solutions for a given value 
of L, only the one minimizing P is selected, which fixes 
the value of n; see below. 

The parameter k is related to the confinement A by 
using Eq. (20): 

A = 2L- [£{ip + £)K, k) + S({n — e)K, fc)], (28) 

where we used the expression (25) for q, qso = sKik), 
and the fact that £ is an odd function. When n = 2p+\ 
such that e = 0, we obtain 

A = 2L - ^£{{2p + l)K, k) = 2L- ^(2p -b l)E[k) 

q q 


■“'-S3' 


(29) 


where E{k) is the complete elliptic integral of the second 
kind [50, p.487]. When n = 2p such that e = 1, we have 


A = 2L- - [£((2p + 1)K, k) -b £(2p - 1)K, k)] (30) 

q 




( 31 ) 


Therefore, the relation between k and A is the same for 
symmetric and antisymmetric solutions. The decrease of 
the applied pressure P as a function of the confinement 
A is given by the parametric equation {A{k), P{k)) given 
by Eqs. (27) and (31) and using Eq. (25). 

Flowever, we still have to determine the optimal value 
of n minimizing P. For L = Ntt (P/A = N/2), it can be 
shown that n = N provided k is small enough. Above 
some threshold, k = k*, we have n = N — 1. As A{k) 
is an increasing function of k [Eq. (31)] this means that 
for small confinement we have n = N before reaching 
a threshold, A(fc*), above which the compressed sys¬ 
tem prefers to remove half a wavelength from the pro¬ 
file, n = N — 1. Flowever, if A(fc*) > Ap, where Ap is 
the critical confinement for which the periodic solution 
is unstable against the localized solution, then n = N 
for all the values of confinement where the periodic so¬ 
lution is stable. The threshold k* can be obtained by 
searching for which value of k, the pressure P has the 
same value for n = N and n = N — We found that 
A{k*)/L ~ {k*)'^ = 27r/3P = A/3L, where we used a first 
order expansion for the confinement which is enough for 
this discussion (using the full expression leads to the same 
conclusion). Comparison with Ap/L = }?jI?, which is 
derived in the next section, shows that for L > 3A the 
optimal value of n is always n = N (with L = Ntt). 
Notice that it leads also to fc < 1/3. Erom now on, we 
assume that the length of the sheet is at least as large as 
three times the wrinkle wavelength. As a result, when N 
is odd, the solution is symmetric and when N is even, it 
is antisymmetric. For L = Nn, we thus have n = N and 
the expression for the pressure reads 


P=^K{kfp 

TT^ 


7r2(l-2A:2) 

AK(kY 


(32) 


Let us summarize the scheme for calculating the ex¬ 
act periodic solution. Given L = Ntt, we have n = N. 
Given A/P we find k form Eq. (31). The values ob¬ 
tained for n and k are substituted in Eq. (25) to obtain 
q. These values of k and q are used in Eqs. (24b) and 
(27) to obtain the height profile and the pressure. This 
solution is unique thanks to the monotonic-increasing 
nature of the right-hand-side of Eq. (31). Finally, the 
pressure-displacement relation is obtained parametrically 
from Eqs. (31) and (32); see Fig. 4. Figure 2 shows the 
evolution with increasing confinement of the exact peri¬ 
odic solution for L = hr (P/A = 3.5) comparing it to 
the numerical solution of Eq. (8). For A/A < 0.2825, the 
numerical solution is periodic. This threshold is close 
to the critical wrinkle-to-fold confinement computed in 
Section IV for large sheets, Ap/A = 0.2857. 

Since it is expected that this exact periodic solution 
is unstable against localization for large enough confine¬ 
ment, we give below its expansion for small A. We 
expand Eqs. (31) and (32) in small k to obtain A = 
Lk^ + 0{k^) and P = 2 — 2fc^ -b 0{k‘^). Thus, to leading 
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order in small relative confinement, we have 

P^2(^l-^y (33) 

The amplitude of the profile is A = 2k/q (see Eq. (24b)). 
For L = Ntt, we have 



The corresponding profiles for small confinement are 

h{s) = A cos s, symmetric profile, (35a) 

h{s) = Asms, asymmetric profile. (35b) 

Thus, to leading order in A/L, our exact solution re¬ 
produces the expressions for P, A, and h{s), known for 
the wrinkle state [2, 6, If, 25]. Notice that, because 
we consider L > 3A, and a transition to the localized 
solution occurs for A/L = X^/L^, the maximum rel¬ 
ative compression relevant for the periodic solution is 
A/L < 1/9 ~ 0.11, which is indeed small. Therefore, 
a first-order expansion is quite sufficient for describing 
the wrinkle state. 


IV. APPROXIMATE LOCALIZED SOLUTIONS 
A. Construction of the localized solution 

We have not been able to find exact physical solutions 
for folds in a finite sheet. Exact localized solutions of 
Eq. (8) for a finite sheet do exist; yet, these profiles are 
not solutions of the present physical problem, as they 
have a finite height at the boundaries. We describe them 
in Appendix A for the sake of mathematical interest and 
possible relevance for other systems to be studied in the 
future. 

We resort to a perturbative calculation, extending the 
multiple-scale analysis of Ref. [44] to a finite system. The 
main motivation is to enable an accurate analysis of the 
wrinkle-to-fold transition, as will be presented in Sec. V. 
We therefore assume a very long (yet finite) sheet com¬ 
pared to the wrinkle wavelength, L 3> A. It has been 
established that, as L is made larger, the region of stable 
wrinkles shrinks and their amplitude diminishes [43, 44]. 
Thus, the critical values of the pressure and displacement 
for the wrinkle-to-fold transition, Pp and Ap, can be as¬ 
sumed arbitrarily close to the ones for the flat-to-wrinkle 
transition, Pw = 2 and Aw = 0, as L is taken to be 
arbitrarily large. 

At the transition from a wrinkle state to a localized 
one, we anticipate that the undulations of wavelength 
A will be attenuated over a much larger length scale of 
order L. Therefore, the sheet profile contains two length 
scales, a short one (A) and a long one (P). To obtain 
this profile near the transition, we then substitute in the 


energy, Eqs. (11)-(13), a multiple-scale function of the 
form 

OO 

h{s) = cos{kcs)''^e^Hj{es) ~ e cos(fccs)Pi(es), (36) 
i=i 

and consider hinged boundary conditions, 

h{±L/2) = h{±L/2) = 0. (37) 

In this expansion, the small parameter e is given by 

e=(Pw-P)'/^ (38) 

where Pw = 2 is the critical flat-to-wrinkle pressure, kf. is 
the wavenumber of the fast oscillations; from the known 
flat-to-wrinkle transition we expect to get kc = 1 [2, 42]. 
We have selected a symmetric profile; an antisymmetric 
one is obtained by replacing the cosine with a sine and 
leads to a similar analysis. For simplicity, we restrict the 
discussion to commensurate sheets, L = Nir {i.e. L/\ = 
N/2), where odd and even N correspond respectively to 
symmetric and antisymmetric solutions. (We will later 
check the effect of this approximation numerically.) 

As in Ref. [44], we substitute the profile Eq. (36) in 
Eq. (13), expand in powers of e, and integrate over the 
fast oscillations (which cancel the terms proportional to 
e and e^), to obtain a systematic expansion of the energy, 

G = Goe^ + Gie* + 0(e®), (39) 

Analysis of the leading order reproduces the known wrin¬ 
kling transition with Pw = 2 and fee = 1 (A = 27r) with¬ 
out any constraint on Pi. This gives Go = 0, reducing 
the energy of the system to G = Gic'^. The function Hi 
is determined by minimizing the functional Gi, given by 

5 IZ Z ^ S"'* + J"') <*> 

where S = es, and a prime denotes a derivative with 
respect to S. Variation of Eq. (40) gives the amplitude 
equation for Hi [44], 

H'/iS) + Ih! - \hi = 0. (41) 

The boundary conditions of vanishing height, h{AL/2) = 
0, are automatically satisfied since cos((A/2)7r) = 0 for 
odd N. The boundary conditions of vanishing bending 
moment, h(±L/2) = 0, impose 

H[{S = ±eL/2) =0. (42) 

Equations (41) and (42) always have the trivial con¬ 
stant solution. Hi = \p2, corresponding to wrinkles. In 
addition, Eq. (41) has solutions in terms of Jacobian ellip¬ 
tic functions [51]. Out of the twelve Jacobian functions, 
only one is found to provide a physical solution, [57] 

Hi{S) = 4:Kdn{KS,m), k =— . (43) 

2-v/2 — m 
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FIG. 3: Localized height profiles for a sheet of length L/A = 
10.5 for various values of the confinement parameter. 


In this section using the other definition of the modulus 
m [51] makes the presentation more elegant. It is related 
to the modulus k of Sec. Ill by m = The modu¬ 
lus m (0 < TO < 1) changes with the displacement A 
(see Eq. (50) below). The pressure is also related to the 
modulus TO using the boundary condition (42),[58] 

neL = 2K{m), (44) 

where K{m) is the complete elliptic integral of the first 
kind [51], which is half the period of the function dn [51], 
k(to) is defined in Eq. (43), and e = y/2 — P. Equations 
(36), (43) and (44) yield the height profile, /i(s), of the 
sheet for a given to. 


h{s) = 4Kedn(Kes, to) cos(s), 
8K{m) ^ f 2K{m) 


V L 


cos(s). 


(45a) 

(45b) 


for L/A = N/2. 

The resulting profiles are demonstrated in Eig. 3. In 
the limit to —>■ 0, the function dn(M,TO) —>■ 1, and 
from Eq. (45b) we recover the wrinkled profile, h{s) = 
V2e cos(s) [43]. In the opposite limit, to —)• 1, we 
have dn(u, to) —>■ l/cosh(u), which recovers the ap¬ 
proximate localized fold for an infinite sheet [43, 44], 
h{s) = 2ecos(s)/cosh(es/2). These two limits are further 
treated in Sec. IV B. Thus, to serves as the order param¬ 
eter of the wrinkle-to-fold transition, as will be discussed 
in Sec. V. The antisymmetric counterpart of Eq. (45b) is 
obtained by replacing cos(s) with sin(s). 

To further characterize the extent of localization, we 
define the decay length of the envelope as. 


gi(0) _ L 

Hi{0) 2y/mK{m)' 


(46) 


where the last equality follows from Eq. (45b). The de¬ 
cay length diverges (^ —>■ oo) when in ^ Q, and van¬ 
ishes (^ —)■ 0) as TO —>■ 1. This defines the two limits 
of weakly and strongly localized profiles to be discussed 


FIG. 4: (Color online) Pressure as a function of displacement 
for a sheet of length L/A = 10.5, showing the second-order 
transition between the wrinkle state (exact solution, green 
solid line) and the fold state (multiple scale approximation, 
red dash-dotted line). P is the physical dimensional pressure. 
Numerical solution (circles) of Eq. (8) indicates that the para¬ 
metric equations (31) and (32) for the wrinkle state, and (47) 
and (50) for the fold state, provide the correct evolution of 
P as a function of A. The pressure-displacement relation for 
a confined infinite sheet, Eq. (58), is shown for comparison 
(exact solution, black dashed line). 


below. Note, however, that the leading order of the 
multiple-scale analysis is not strictly valid in the strongly 
localized limit. The crossover between these two limits, 
^ = L, where the localization should become observable, 
occurs at TO ~ 0.096. 


B. Pressure, displacement, and amplitude 

Using Eq. (44) together with Eq. (38), we obtain the 
expression for the pressure associated with the localized 
solution, 

P(to) = 2 — ^(2 — TO)iL(TO)^. (47) 

Since we know that for an infinite sheet the wrinkle state 
is always unstable against the localized state [44], we ex¬ 
pect that this transition should occur at an arbitrarily 
small displacement for a sufficiently long sheet. There¬ 
fore, we expand the expression (12) of the displacement 
to the leading order in e, 

A ~ / dsh'^. (48) 

2 J-L/2 

Using the expression (45a) for /i, averaging over the fast 
oscillations (short length scale), and performing the inte¬ 
gral over the slow envelope (long length scale), we obtain 

A = 8Ke£’(KeL/2, to), (49) 

where £(x,m) is the Jacobi epsilon function, and k is 
given as a function of to in Eq. (43). Using Eq. (44), 
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the displacement associated with the localized solution 
reads, 

16 

A(m) = 8KeS{K{m),m) = —K[m)E{m), (50) 

1j 

where E{m) is the complete elliptic integral of the sec¬ 
ond kind [51]. Since both A, given by Eq. (50), and P, 
given by Eq. (47), are functions of m, the evolution of 
the pressure as a function of the displacement is given 
by the parametric form (A(to), P{m)). This evolution is 
demonstrated in Fig. 4. 

1. Weakly localized limit 


used to experimentally characterize the evolution of the 
localized profile [42]; see also Sec. VI and Fig. 7. (The lo¬ 
calization makes the second extremum shift slightly from 
s = TT, but this shift is of higher order than our present 
approximation.) From Eqs. (45b) and (52) we get 

, (56a) 

^ + 2V2 cos ypj , (56b) 

which are valid for A > Ap. The slopes of both Aq and 
Ai diverge at Ap, with an exponent of 1/2, as is typical 
for a mean-field second-order transition. 


As m grows from zero the profile ceases to be purely 
periodic and begins to localize. In the weakly localized 
limit, m ^ 1, the symmetry has already been broken, 
but the localization length is still larger than the system 
size, ^ > L. In this limit, therefore, the deviation from 
the wrinkle state will be hard to resolve experimentally. 
Nevertheless, as mentioned above, m ~ 0.1 is sufficient 
to reach L. 

The displacement given by Eq. (50) is an increasing 
function of m. We notice that, in the limit m —>■ 0, the 
displacement takes a finite value, 

Ap = A(0) = 47rVA = \^/L. (51) 


This means that the localized solution emerges only be¬ 
yond this finite displacement in the case of sheets of finite 
length. For a confinement smaller than Ap, the shape 
of the sheet is described by the periodic solution con¬ 
structed in Sec. III. In Appendix C we obtain the critical 
wrinkle-to-fold confinement, up to a constant prefactor, 
Ap ~ A^/L, from a simple scaling analysis. 

For small m we have 

2 2 
TT^ID ^ 

A{m) = Ap ++ 0{m^), (52) 

P(m) = 2 - ^ + Oim^). (53) 

The parameter m can be eliminated to obtain 


P~2- 


2Ap 

~L~ 


-6 


(^) 


(54) 


The profile for m = 0 is obtained from Eq. (45b), 

47r 

h(s) = —coss, TO = 0. (55) 

Jj 

Although this profile describes periodic wrinkles, it has a 
finite amplitude since it is obtained for a finite displace¬ 
ment, A = Ap. 

Once TO > 0, the profile is no longer periodic and its 
extrema have different heights. The height of the cen¬ 
tral fold, Af) = |h(0)|, and that of the extremum next 
to it, Ai ~ |h(7r)|, are of special interest since they are 


2. Strongly localized limit 


In this limit of to 1 we have 

A(to) ~ ^K{m), (57a) 

Jj 

1 ft 

P{m) ~2 —-^K{mY. (57b) 

Eliminating the parameter to between these two last re¬ 
lations leads to 


P = 2- 


16 




(58) 


which coincides with the exact pressure-displacement re¬ 
lation for infinite sheets [45]. In this limit, we also have 
dn(M, to) —)■ l/cosh(u) [51], such that the profile (45b) 
becomes 


h{s) 


A cos(s) 

2 cosh(s/^)’ 


e = 8/A, 


(59) 


which coincides with the approximate localized fold for 
an infinite sheet [43, 44]. The amplitude Aq is thus given 
by the following simple relation. 


Ao = A/2, 


(60) 


which again is identical to the exact result for an infinite 
sheet [45]. 


V. WRINKLE-TO-FOLD TRANSITION 

Below a critical value of eL = epL = 2^/^7r, Eq. (44) 
has no solution for to. Thus, below this critical confine¬ 
ment, the localized solution does not exist, and the only 
possible envelope is the trivial constant one, correspond¬ 
ing to wrinkles. At epL itself Eq. (44) admits to = 0 as a 
unique solution, and for larger values of eL, with to > 0 
as a unique solution, the envelope grows continuously, 
corresponding to increasingly localized patterns. Accord¬ 
ingly, we define the order parameter of the wrinkle-to-fold 
transition as to, and its control parameter as 

r = eL, Tp = epL = 2^/^7r. 


(61) 










FIG. 5: (Color online) Phase diagrams of an incompressible fluid-supported sheet, (a) Phase diagram on the displacement- 
length plane. For 0 < A < A^/L the sheet forms extended wrinkles, and for A > A^/L a localized fold is formed. The 
rectangle shows the range of the experimentally determined phase boundary [42] (see also Fig. 7). (b) Phase diagram on the 
pressure-length plane. P is the physical dimensional pressure. Wrinkles are stable for 2 — 2 \^ fL^ < P < 2 , and a fold is stable 
for P <2 — 2')? j I? . The wrinkle-fold line is a second-order transition. Note how the region of stable wrinkles in both diagrams 
vanishes in the limit of an infinite sheet. In both diagrams, the theoretical transition curves are compared to the numerical 
transition obtained by solving Eq. (8) and analyzing the periodicity of the numerical solution from the relative height of its 
extrema. Numerical results for both commensurate [L being an integer or half-integer multiple of A) and incommensurate 
sheets are presented, showing a negligible discrepancy between these two and the theoretical prediction for L > 2A. 


A. Landau expansion 


We repeat here these central results: 


A clear way to present the continuous transition is 
through a Landau expansion of the energy, G, in small 
order parameter. The Landau expansion is to be per¬ 
formed on the energy prior to minimization. For this 
purpose, we use the following form for the envelope: 

HI {S) = ^ dn S, m'] , (62) 

V2 — m \ T ) 

where we have substituted the boundary condition (44) 
in the argument of the profile given by Eq. (43), but not 
in its amplitude. Therefore, H*{S) satisfies the boundary 
conditions but does not yet minimize the energy. 

Substituting Eq. (62) in Eq. (40) and expanding to 
fourth order in small m, yields 

AG ~ (r| - t^)(to^ -b m^) + 7^(35r| - (63) 

iZo 

where AG = 32L(G — Gw)/e^, and Gw = Le'^/4 is the 
energy of the wrinkles. The appearance of the term 
has no significance for the transition because of its van¬ 
ishing prefactor. Upon minimization of AG with respect 
to m we obtain, 


m ~ 8 [(r - rF)/3rF]^^^ , (64) 

with a critical exponent, /3 = 1/2, as usual for a mean- 
field second-order transition. 

Using the expression (61) of rp in Eqs. (33) and (38) 
for P and A, we recover the critical values derived in 
Sec. IV B1 from the analysis of weakly localized folds. 


Pp = Pw-r|/L2 =2-87rVL^ 

Af = t|/(2L) = (66) 

The value of Ap — 39.5/L corrects an earlier, higher es¬ 
timate of Ap = A8jL, which was obtained by using an 
Ansatz profile that was not an energy minimizer [43]. 
Equations (65) and (66) confirm that the wrinkle-to- 
fold transition can be made arbitrarily close to the flat- 
to-wrinkle one with increasing size of the sheet. Con¬ 
sequently, this separate second-order transition appears 
only in finite systems. 

In Fig. 5 we summarize the phase diagrams of the in¬ 
compressible fluid-supported sheet on the A-L and P~L 
planes. Note that these scaled diagrams are parameter- 
free. 


B. Matching at the transition 


Below the wrinkle-to-fold transition and considering 
small wrinkles (to first order in A/L), we have found in 
Sec. me, 


P = 2 




h{s) = Aq cos(s). 


(67) 

At the wrinkle-to-fold-transition, we have shown that 
A = Ap = AVL and P = Pp = 2 - 2X‘^/L‘^. Slightly 
above the transition, the weakly localized solution is 
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- Wrinkles -Fold (i = 5.5A) - Wrinkles - AJ\(L = 7.5\) - AJ\(L = ^.5X) 

Fold {i= 10.5A) -- Fold(L = 3.5A) A„/A (i = 10.5A) AJX[L = 7.b\) -- A„/A{i = 3.5A) 

- Fold{i = 7.5A) ® Numerics . Aj/A (i = 10.5A) - A„/A (i = 5.5A) - Aj/A{i = 3.5A) 



FIG. 6: (Color online) (a) Evolntion of the rescaled pressnre as a fnnction of the relative displacement for various sheet lengths 
L. P is the physical dimensional pressure. The graphs show the evolution of the pressure for the periodic state (wrinkles) 
together with the bifurcations toward the localized state (fold) at A/L = Ap/L = jl? and P = 2 — 2Af/L. These 
bifnrcations decrease the energy of the system compared to the periodic state, (b) Evolution of the amplitndes of the pattern, 
rescaled by the wrinkles wavelength, as a function of the relative displacement for the same sheet lengths as those used in panel 
(a). The wrinkles amplitude grows as until the points A/L — }? jlP and A/A = XjLn is reached, where the localized 

state emerges. The evolution of the amplitude of the central fold, Aq, and of the fold next to it, Ai, are shown. 


characterized by (see Sec. IVB 1), 


P = 2 




- 6 


A-Af\ 


Ao = ^+2V2 

1j 

h{s) = ^ cos(s), 

1j 


at the transition 


The measured quantity is the height profile, including 
the wrinkles wavelength. A, and height extrema [42]. 

The order parameter m, which we have defined here, 
characterizes the decay of the envelope away from the 
center of the fold. Experimentally, one may look instead 
at the height difference between two consecutive local 
extrema Aq and Ai [42]. These two definitions of the 
order parameter are equivalent close to the transition. 
Indeed, from Eqs. (52) and (56) we find 


Comparing Eqs. (67) and (68), we see that the tran¬ 
sition is continuous. The discontinuity appears in the 
pressure derivative, with {dP/dA)^^^- = —2/L and 
{dP/dA)^^^+ = —djL. Hence, the continuous transi¬ 
tion is of second order (discontinuity in the second deriva¬ 
tive of E with respect to A). The pressure associated 
with the localized solution is smaller than the one asso¬ 
ciated with the periodic solution once A > Ap, as illus¬ 
trated in Fig 6a. Since P = dE/dA, the localized state 
has a smaller energy compared to the energy of the peri¬ 
odic state. The latter is thus unstable once A > Ap. The 
evolution of the amplitudes of the wrinkle and localized 
states as the displacement increases is shown in Fig 6b. 


VI. DISCUSSION 
A. Experimental consequences 

The system parameters are the length L and bending 
modulus B of the thin elastic sheet, and the density p 
of the liquid. The experimental control parameter is ei¬ 
ther the displacement A or the force per unit length P. 


A7exp = 1 - 4^ = ^[sin(A7r/2L)]^TO -k 0{rrP), 
Aq Z 


8L2 


. + 0{m^,X^/L^). 


(69) 


Using Eq. (52), or equivalently Eq. (64), this result is 
written as a function of the experimental control param¬ 
eter A: 


TT^ A /A-Ap 

•\/2 L \ 


1/2 


(70) 


where Ap is given by Eq. (51). 

The evolution of the pattern amplitude as a function 
of the displacement can also be measured. Such data are 
available in Ref. [42], whose main purpose was to show 
the universality of the folding state where all the data col¬ 
lapsed onto a master curve for large enough confinement. 
Indeed, in the large-confinement regime, the sheet defor¬ 
mation is localized in a small region, comparable to the 
wrinkle wavelength A, and does not depend significantly 
on the sheet length L, provided that L is sufficiently large 
compared to A. The data of Ref. [42] are reproduced in 
Fig. 7 together with the evolution of the amplitudes pre¬ 
dicted by the present theory. (Note that, to comply with 
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the presentation of the experimental data, the amplitude 
in Fig. 7 is drawn as a function of A/A rather than A/L 
as in Fig. 6b.) 

Our theory is valid for large enough sheets and con¬ 
sequently for a sufficiently small relative confinement 
A/L ^ 1. See, for example, the deviation of the ampli¬ 
tude Ai at large A/L for a small system with LJX = 3.5 
(Fig. 7). On the other hand, at large A/A, as the lo¬ 
calization length becomes significantly smaller than L, 
the exact predictions for an infinite sheet [45] become 
accurate for finite sheets as well (see Fig. 7). Finally, 
for the localized state we have assumed a commensurate 
sheet {L being an integer or half-integer multiple of A). 
Yet, as shown in Fig. 5, the effect of incommensurabil¬ 
ity becomes negligible already for sheets larger than 2A, 
making it inconsequential experimentally. 

The authors of Ref. [42] inferred from their experi¬ 
ments that the wrinkle-to-fold transition occurred at a 
value of A ~ 0.3A for all sheet lengths L, whereas our 
theory gives an L-dependent critical displacement (see 
Fig. 5a). As seen in Fig. 7, it was natural to draw that 
conclusion given the experimental error. The rectangle 
shown in the diagram of Fig. 5a, representing the range 
of the experimentally determined thresholds (Fig. 7), is 
consistent with the theoretical prediction. Clearly, accu¬ 
rate experimental data focusing on the transition region 
are still needed. For this purpose we summarize the pre¬ 
dictions for the amplitude evolution below and above the 
transition: 


For A < Ap : 

A 1 /aV^^ 


A TT \l) ’ 

For A > Ap : 

Ag _ A / A — A] 

A ttL t!" V L 

Ai X y/2 /27r2 


1/2 


A-Ap 


1/2 


(71a) 

(71b) 

(71c) 


In Sec. IV, we derived the properties of the critical 
wrinkle-to-fold transition at which the localization length 
^ diverges. As for critical phenomena in a finite system, 
the transition will be observable in practice only when the 
localization length becomes smaller than the system size, 
^ < L. We find, however, that the displacement required 
to get ^ ~ L is larger than Ap by only 0.4%. Thus, the 
finite size of the system should have a negligible effect on 
the experimentally observed critical displacement. 



FIG. 7: (Color online) (a) Definitions of the amplitudes Ag 
and Ai for the localized state, (b) Comparison between the 
experimental evolution of Ag and Ai (rescaled by A) with in¬ 
creasing confinement for finite sheets [42], and the prediction 
of Eq. (45b) where Ag = |h(0)| and Ai = |h(7r)|. The evo¬ 
lution of the amplitudes of the exact solution for an infinite 
sheet is also presented [45]. 


A higher-order theory (next order in e of the multiple- 
scale expansion) would not change the main conclusions 
obtained above. It is expected to slightly modify some 
quantitative predictions such as the value of the critical 
wrinkle-to-fold confinement, Ap. The possibility to dis¬ 
criminate experimentally between the present theory of 
the transition and a higher-order one is doubtful. 

Although we found exact localized solutions for the 
governing equation (8), these do not give the physical 
profile of the sheet (see Appendix A). Mathematically, 
the inability to derive physical localized profiles from 
these solutions may indicate that such profiles belong to 
a different branch of solutions to the non-linear governing 
equation. Such a branch is yet to be found. If it exists, 
it will probably coincide, in the limit L —>■ oo, with the 
exact solution known for this limit [45]. 

An important extension to the present theory is the 
inclusion of finite compressibility of the sheet. This will 
allow us to treat the flat-to-wrinkle transition, and thus 
construct the complete “phase-diagram” of the floating 
sheet including the two transitions. This extension will 
be presented in a forthcoming publication. 


B. Conclusion 

In the present work we have derived an exact solu¬ 
tion for the wrinkle state in a finite fluid-supported sheet. 
The theory developed for fold localization, however, is re¬ 
stricted to large sheets; it is correct to first order in X/L. 
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Appendix A: Exact localized solutions to Eq. (8) 


In this Appendix we derive exact localized solutions 
to Eq. (8), which are not physical solutions of our 
problem treated here, but may be of use in future re¬ 
search. Equation (8) is a member of the sine-Gordon- 
modified-Korteweg-de Vries hierarchy [52] in which the 
sine-Gordon equation is a lower member. Consequently, 
a solution of the sine-Gordon equation is also a solution 
of Eq. (8). Using the separation of variable proposed by 
Lamb [53], the following function 

= 4arctan(F(s)/G(t)), (Al) 



is a solution of the sine-Gordon equation provided F and 
G satisfy the following differential equations [54] 

(F')^ = -kF^ + ^lF‘^ + A, (A2a) 

{G'f = kG^ + 1)G2 - A, (A2b) 

where k, fj, and A are arbitrary constants. The function 
(Al) with t = s is also a solution of Eq. (8) if F and G 
satisfy the relations (A2) and provided that 


FIG. 8: (Color online) Evolution of the normalized height of 
the exact solution (A6a) at the boundaries as a function of the 
normalized displacement. The dots indicate the position of 
the threshold, Ap/A, below which the profile is periodic. The 
shaded area corresponds thus to the wrinkle state. This area 
is delimited by the relation h{±L/2)/h{0) = 0.279 Ap/A = 
0.279 A/L. 


each parity of </>) as a function of A and (3 as follows: 


P = 2 - 4/r. (A3) 

Introducing some scaling parameters [55], Af{j5s) = F{s) 
and g{ujs) = 1/G(s), Eqs. (A2) become 

U'f = [-nA^f^ + + AA-2] , (A4a) 

{g')^ = w ^ [—Xg'^ + (m ~ 1 ) 5 ^ + ■ (A4b) 

Solving Eqs. (A4) allows to determine / and g and to 
obtain an exact solution of Eq. (8). Those equations are 
satisfied by the Jacobi elliptic functions with arguments 
/3s and ujs and parameters kf and kg. The various pa¬ 
rameters are fixed by comparison between Eqs. (A4) and 
those satisfied by the Jacobi functions: 

sn(s, k) : (j/')^ = k'^y'^ — (1 + + 1 (A5a) 

cn(s, k) : {y'f = -k^y'^ + (2fc^ - l)y'^ + l-k^ (A5b) 

dn(s, k) : {y'Y = —y^ + (2 — k'^)y^ — l + k"^ (A5c) 


Hinged sheets 

We choose / = cn(/3s,A:/) and g = sn(a;+s, fc+) or 
g = cn(a;“s, k~) such that 

(j){s) = 4arctan[Acn(^s, A:/) sn(a;'''s, fc^)], (A6a) 

(j){s) = 4arctan[Acn(/3s, A:/) cn(a;“s, fc“)]. (A6b) 

Equations (A6a) and (A6b) lead to symmetric and an¬ 
tisymmetric solutions respectively. Comparison between 
Eqs. (A4) and Eqs. (A5) for that particular choices of / 
and g fixes six parameters (among eight parameters for 


w+ = 


[1 + (1 + A 2 )/ 32 ] 


1/2 


k+ = A 


kg = A 


1-h A2 

^ +13^ 

1 + A^ ^ 

l-hA2 


- 1 


1/2 


-,1/2 


- 1 


1 + (1+A2)^2 

/32-^^2(/32_ 1) -|l/2 

l-A4-h(l-hA2)2/32 


kf = 


{1 + A'^+p-'^) 


^ 1 + A^ 

2A^ + (A^ - 1)^2 


-2\1/2 


M = 


(1 -h A2)2 


(A7a) 

(A7b) 

(A7c) 

(A7d) 

(A7e) 

(A7f) 


The expressions of k and A are not needed in the fol¬ 
lowing and are not written. The amplitude A of the 
solutions (A6) is related to the displacement A through 
Eq. (2) whereas [3 should be fixed to satisfy the boundary 
conditions. 

For hinged boundary conditions both h and h = (j) cos (/> 
vanish at s = ±L/2. From Eq. (7), this is equivalent to 
(^(±L/2) = ^'(±L/2) = 0. The condition 0(L/2) = 0 for 
the symmetric solution (A6a) leads to the constraint 


dn(a:, A:/) sn(a:, A:/) sn(?/, A:+) y 
cii{x,kf)cn{y,kg)dn{y,kg) x" 


(AS) 


where x = /3T/2 and y = w+L/2. Due to the definite 
parity of the solution, (j>{—L/2) = 0 is automatically sat¬ 
isfied. A similar relation is easily obtained for the an¬ 
tisymmetric solution (A6b) which is not written here. 
For a given amplitude. A, and a given sheet length, L, 
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FIG. 9: (Color online) Evolution of the solution (A6a) for 
different values of A/A and L/A = 5.5, where a:(s), h{s) are 
given by Eqs. (1) and </(s) is given by Eq. (A6a). The cor¬ 
responding evolution for the numerical solution of Eq. (8) is 
also shown. For A/A < Ap/A = 0.182, the numerical solu¬ 
tion is periodic and coincides with the solution reported in 
Section III. 


Eq. (A8) fixes /3. The displacement, A, and the pressure, 
P, are then computed from Eqs. (2) and (A3) respec¬ 
tively. Finally, P is computed from Eq. (9). 

It is however impossible to satisfy both boundary con¬ 
ditions with the solutions (A6). For the exact periodic 
solution (18), both and (j){s) are proportional to 
cn(( 7 (s -|- So), A:) such that they can simultaneously van¬ 
ish at the boundaries with a suitable choice of q (and 
So), see Eqs. (22). Here, ^{L/2) and (j){L/2) cannot 
simultaneously vanish for the same value of /?. There¬ 
fore, the height of the profile assumes a finite value at 
the boundaries as shown in Fig. 8. Notice however that 
/i(±L/2)//i(0) decreases as L increases. 

Nevertheless, despite the fact that it does not properly 
satisfy the boundary conditions, this solution is a very 
good approximation of the numerical solution satisfying 
both boundary conditions. Figure 9 shows a comparison 
between this exact localized solution and the numerical 
solution of Eq. (8). The agreement is good especially 
near the central fold and for large enough confinement. 


Appendix B: Accuracy of Ap and Pp 

The critical displacement, Ap, and pressure, Pp, at 
which the transition from wrinkles to fold occurs has 
been computed in Sections IV B1 and V. These expres¬ 
sions are obtained from matching, at the first order in 
A/L, between the exact periodic and the approximate 
localized solutions, see Section VB, and converge toward 
the exact values in the limit P —>■ oo. Although a rigor¬ 
ous derivation of the error on these two quantities is only 
possible if their exact expressions are known, an estima¬ 
tion of the rate of convergence can be obtained from the 
exact periodic solution. 

Indeed, since at the transition P(A) coincides for both 


FIG. 10: (Golor online) Evolution of 5 a and 5p as a function 
of the sheet length (normalized by the wavelength) obtained 
by solving numerically Eqs. (Bl). The asymptotic behaviors 
(B2) are also shown. 


the periodic and the localized solutions, the exact ex¬ 
pressions of Ap and Pp should satisfy Eqs. (31) and (32) 
for the same k = kc- This is not the case with the ap¬ 
proximate expressions we have obtained in this paper. 
Therefore, we define two quantities, 6a and dp, as fol¬ 
lows 


6a = with P{kc) = Pp, (Bla) 

A (fee) 

i5p = — ^ ^ with A(fcc) = Ap, (Bib) 

P(fcc) 1 c; F, I ; 


where A(fc), P{k), Pp and Ap are given by Eqs. (31), 
(32), (65) and (66) respectively. With the exact expres¬ 
sions of Ap and Pp, we have kc = kc and i5a = <5p = 0. 
The quantities 6a and 5p are thus a measure of the error 
introduced and of the rate of convergence toward the ex¬ 
act expressions of the critical displacement and pressure. 

Expanding both A(fc) and P{k) up to k"^, we obtain at 
the leading order in A/P: 




A^i 

P4 


u r 3A4 
and 5p = —-bO 


A® 

P6 


■ (B2) 


Figure 10 shows that the evolution of 6a and dp as a func¬ 
tion of P/A obtained by solving numerically Eqs. (Bl) 
agrees well with the asymptotic expression (B2) even 
when P/A is not so large. Therefore, the expressions 
of Ap and Pp derived in the paper converges toward the 
exact values as L~^ and P“^ respectively. 


Appendix C: Scaling approach 

In this section, we show that some of the results ob¬ 
tained in the main text, and in particular the critical 
confinement at which a wrinkle-to-fold transition occurs, 
can be qualitatively recovered using a simple scaling ap¬ 
proach. As formulated in Sec. II, the bending energy of 
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Wrinkles Fold 



FIG. 11: Schematics for wrinkles and fold patterns. 


the sheet E\j and the deformation energy of the substrate 
Eg read 



WB 


(Cl) 

Eb = 

~ o i 

1 jXds, 


2 J 

-L/2 



Wpg 

nL/2 

/ hJ cos (j>ds, 

(C2) 


2 ^ 

'-L/2 



and the displacement along the direction of confinement 
is given by 

.L/2 

A = / (I — cos 4>)ds. (C3) 

J-L/2 

1. Small displacement 

For small displacement, leading to small sheet defor¬ 
mations 1), wrinkles emerge with an amplitude A, 
wavelength A and curvature 0 ^ A/X^, see Fig. 11. The 
energies thus scale as 

Eb - BWLA^/X'^; E^ ~ pgWLA^. (C4) 

The balance of these two energies leads to 

A ^ {Blpgf/X (C5) 

In the limit ^ <g; 1, the relation between the amplitude 
and the displacement is given by 

rL/2 t-L/2 

A~ / {dh/ds)"^ ds ^ L{A/Xf . (C6) 

J-L/2 J-L/2 


2. Large displacement 

The total energy of the system can be written as 


f ; 

.L/2 ^ 

jJ^ds + 

P_9 

hJds 

(C7) 

2 7 

-L/2 

^ J-L/2 



.L/2 

/ h^{l- 

COS J) )ds 


(C8) 

2 J 

-L/2 




RAt/^ -1- pgAA^ 

- pgA^A 

1 

(C9) 

pg 

X'^kjX E KA 

^ - A^A 

1 

(CIO) 


where A is the spatial extent of the localized pattern, 0 
its typical curvature and A its amplitude. 

For wrinkles, we have K ^ (j) ^ AjXJ and A ^ 

see Eq. (C6). The wrinkles energy can thus be 
written as 

Ew/W = pgX^ [aA/X-b{X/L){A/Xf] , (Cll) 

where a and b are numerical constants. For a fold, we 
have A ~ the fold width, cj) ^ £~^, and A ~ A, see 
Fig. 11. The fold energy can thus be written as 

Ef/W = pg [(cV4c)AV^ + c£A 2 - dA^] , (C12) 

where c, c and d are numerical constants. A minimization 
of Ep with respect to i gives £ = (c/2c)A^/A and the 
energy becomes 

Ef/W = pgX^ [c A/A - d(A/A)3] . (C13) 

Requiring that the wrinkles and fold energies match at 
A = Ap gives a quadratic equation for Ap. Demand¬ 
ing that Ap vanishes when L diverges yields a = c. 
A wrinkle-to-fold transition occurs when Eyp > Ep or 
equivalently, using a = c, when A > Ap = {b/d)X^/L. 
Thus, this simple analysis has recovered the correct scal¬ 
ing of Ap with respect to A and L. As shown in Sec. V, 
the prefactor b/d turns out to be exactly 1. 
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